Inversion of DNA charge by a positive polymer via fractionalization of the polymer 

charge 



T. T. Nguyen and B. I. Shklovskii 

Theoretical Physics Institute, University of Minnesota, 
116 Church St. Southeast, Minneapolis, Minnesota 55455 

Charge inversion of a DNA double helix by an oppositely charged flexible polyelectrolyte (PE) 
is widely used for gene delivery. It is considered here in terms of discrete charges of DNA. We 
concentrate on the worst scenario case when in the neutral state of the DNA-PE complex, each of 
the DNA charges is locally compensated by a PE charge and show that charge inversion exists even 
in this case. When an additional PE molecule is adsorbed by DNA, its charge gets fractionalized 
into monomer charges of defects (tails and arches) on the background of the perfectly neutralized 
DNA. These charges spread all over the DNA eliminating the self-energy of PE. Fractionalization 
leads to a substantial charge inversion of DNA. We show that fractionalization mechanism charge 
inversion works also for nonlinear polymers such as dendrimers. Remarkably, such fractionalization 
happens for adsorption of both PE or dendrimers on a two-dimensional charged lattice, as well. 
Relation of fractionalization to other mechanisms of charge inversion is discussed. 

PACS numbers: 87.14.Gg, 87.15. Nn, 87.16.Dg 



I. INTRODUCTION 

Inversion of the negative charge of a DNA double helix 
by its complexation with a positive polyelectrolyte (PE) 
is used for the gene delivery. The positive charge of DNA- 
PE complex facilitates DNA contact with a typically neg- 
ative cell membrane making penetration into the cell hun- 
dreds times more likelyEl. Charge inversion of DNA-PE 
complexes was confirmed recently by electrophoresis!! If, 
at a given concentration of long DNA helices, the concen- 
tration of shorter PE molecules increases, at some criti- 
cal point the electrophoretic mobility of a DNA-PE com- 
plex changes sign from negative to positive. Intuitively, 
one can think that when a PE completely neutralizes a 
DNA double helix new molecules of PE stop adsorbing 
on DNA. Indeed, the Poisson-Boltzmann approximation 
for description of screening of a DNA helix by any coun- 
tcrions including PE does not lead to charge inversion. 
The counterintuitive phenomenon of charge inversion of a 
macroion by oaap^te^ch.ar g ed r E' [ E 1 has attracted signifi- 
cant attentio Jflra#M»rjy. It can be explained 
if one takes into account that the surface potential of an 
already neutralized DNA is locally affected by a new ap- 
proaching PE molecule, or in other words —taking into 
account correlations between PE moleculesOO. Due to 
repulsive interaction between PE molecules a new PE 
molecule pushes aside PE molecules which are already 
adsorbed on DNA surface and creates on the surface an 
oppositely charged image of itself. The image attracts 
the new PE molecule leading to charge inversion. This 
phenomenon is similar to attraction of a charge to a neu- 
tral metal. 

For quantitative consideration, charges of DNA are 
often assumed to be smeared and to form uniformly 
charged cylinder. This approach seems to be justified 
when density of charge of PE is larger than density of 



charge at the DNA surface so that most of DNA surface 
is empty. However, it is clearly far from satisfactory when 
these densities are almost equal and PE charges strongly 
compete for charges of DNA (see figures below). Ap- 
proximation of uniform charge also ignores interference 
between chemical structure of DNA surface and of PE. 
Therefore, generally speaking, it is not even clear whether 
charge inversion exists in the case of discrete charges or it 
is just an artifact of the assumption of uniformly smeared 
charge. In this paper, we consider effects of discreteness 
of — e charges of DNA. We show that in this case charge 
inversion exists as well. It can be explained as a result of 
the "fractionalization" of charge of PE molecules. Such 
explanation turns out to be even simpler and more visual 
than for the model of smeared charges of DNA. 

Negative elementary charges of DNA phosphates are 
situated along the two spirals at the exterior of both he- 
lices. When unfolded, each spiral is an one-dimensional 
lattice of such charges, with the lattice constant a=6.7A. 
Let us consider a toy model of a PE as a freely jointed 
chain of Z small +e monomers. To maximize the role of 
discreteness of DNA charge we begin from the assump- 
tion that the PE bond length b is exactly equal to the 
distance a between negative charges of a spiral. We call 
such PE "matching" . We also assume that minimal dis- 
tance, d, between a PE charge and a charge of DNA is 
smaller than a. Then PE molecules can attach to a DNA 
charge spiral in such a way that every charge of the spi- 
ral is locally compensated by a PE charge and, therefore, 
DNA is completely neutralized. 

The case of a very short polymer (oligomer) with Z = 3 
is shown in Fig. [I] simplest illustration. The neu- 

tralization by a matching PE is so perfect that it is diffi- 
cult to imagine how another PE molecule can be attached 
to DNA. Thus, it seems to be impossible to overcharge 
DNA. In this paper ,we show that even in this worst pos- 
sible for charge inversion scenario, there is a mechanism 
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FIG. 1: The origin of charge fractionalization. a) One of the 
spirals of negative charges of DNA (white spheres) is com- 
pletely neutralized by positive PE molecules with Z = 3 Their 
charges are shown by black spheres, neutral atoms are not 
shown. A new PE molecule is approaching DNA. b) The new 
PE molecule is "digested", its charge is split in +e charges of 
two tails and an arch (center). 



which brings an additional PE to the neutralized DNA 
and leads to charge inversion. We call this mechanism 
fractionalization and Fig. [l] shows how this mechanism 
works for the case of Z = 3 When a new PE comes to the 
DNA double helix which is already neutralized by PE, it 
creates a place for itself or, in other words, the oppositely 
charged image in the following way. 

In each of Z already-adsorbed PE molecules one PE 
monomer detaches from DNA surface. This leads to for- 
mation of positive defects (tails and arches) and Z nega- 
tive vacancies on DNA. All Z vacancies can join together 
and form a large vacancy of a length Z by shifting of 
adsorbed PE molecules along DNA. A new PE molecule 
is accommodated in this vacancy. As a result of con- 
sumption of this molecule Z defects with charge +e each 
appear on top of the completely neutralized spiral (see 
Fig. 0b). 

This effectively looks as cutting of the new PE molecule 
into Z individual monomers and spreading them out 
along the spiral. In other words, charge inversion of 
DNA happens by fractionalization of the PE molecule 
charge. Of course, none of the chemical bonds is really 
cut, and this phenomenon is solely due to the correlated 
distribution of PE molecules, which avoid each other at 
the DNA spiral. In this sense, fractionalization we are 
talking about is similar to what happens in Pactional 
quantum Hall effeclJi3 or in the polyacetylenetj, where 
many-electron correlations result in the fractionalization 
of the electron charge. 

Fractionalization is driven by elimination of the self- 
energy of free PE molecules. By the self-energy we mean 
the energy of repulsive interactions of Z positive charges 
of the PE molecule in extended conformation which it 
has in the solution. In the fractionalized state, charges 
of monomers are very far from each other and practically 
do not interact, so that the positive PE self-energy is 
eliminated and, therefore, gained. 

In the next section we calculate fractionalization in- 
duced charge inversion by a matching flexible PE. In Sec. 



Ill we discuss what happens when PE does not match the 
DNA spiral of charges so that linear densities of charge 
are different. In Sec. IV we generalize these ideas to ad- 
sorption on two-dimensional lattices of discrete macroion 
charge. We show that in the case of matching flexible PE 
fractionalization works perfectly even in two-dimensions. 
This is interesting because many other physical exam- 
ples of charge fractionalization do not work beyond one 
dimension. Furthermore, this is the first classical ex- 
ample of real two-dimensional fractionalization. In Sec. 
V we generalize our theory to flexible polymers, which 
do not have linear structure. We concentrate there on 
charge inversion of DNA by dendrimers and show that 
in this case fractionalization ideas lead to charge inver- 
sion, too. In Sec VI, we discuss additional mechanism 
of charge inversion related to the fact that DNA charges 
can be accessible from two opposite sides. We conclude 
in Sec. YII. A short version of this paper is published 
elsewhereEJ. 



II. FRACTIONALIZATION INDUCED CHARGE 
INVERSION: A MATCHING 
POLYELECTROLYTE 

Let us now calculate the linear density of the net charge 
of DNA, rj* , using the fractionalization mechanism. The 
chemical potential of the PE absorbed at the spiral is 



Ms = Zk B T\n{rf /r, ) + Zeif>(0) 



(1) 



The first term in the right hand side of Eq. ([!]) is the 
chemical potential of the one-dimensional gas of defects 
(—770 — 0.6e/A is the bare charge density of DNA). We 
used expression for the chemical potential of an ideal gas 
because the Coulomb interaction energy between defects 
at the a distance of a few a is much smaller than fc^T (a ~ 
Ib, where l B — e 2 / DksT ~ 7A is the Bjerrum length.) 
The second term in the right hand side of Eq. (|l[) is the 
repulsion energy of the new PE from the inverted charge 
of the DNA. In this term, tp(0) is the averaged surface 
potential of the DNA helix. We assume in this paper that 
the net charge of DNA is screened by a monovalent salt 
at the screening length r s , which is much larger than a. 
Then ^(0) can be calculated as the surface potential of a 
cylinder with radius of DNA helix R and linear density 
of charge rj* 
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To find rf in the equilibrium state, one has to equate the 
chemical potential of adsorbed PE molecules with that 
of a free PE in the solution. The later one can be calcu- 
lated as following. Due to the repulsive Coulomb inter- 
action between monomers, a free PE in the solution has 
an extended shape to minimize its energy. Therefore, the 
chemical potential of a free PE in solution can be written 
as the self-energy of a rigid rod with the length Za and 
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the linear charge density e/a plus the standard ideal gas 
contribution kgTlD.(Nvo) (N is the number concentra- 
tion of free PE in solution and vq is the volume of a PE 
molecule) : 



Mo = (Ze 2 1 Da) ln(£/ a) + k B Tln(Nv ) , 



(3) 



where C = mm(r s , Za) and D is the dielectric constant 
of water. 

Equating the chemical potentials of Eqs. (||) and ([!]), 
one has 



V(0) = -f In- 
Da a 
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In Eq. (Q) one can interpret the right hand side as a 
"correlation" voltage that (over)charges the DNA to the 
potential "0(0)- Complete analysis of Eq. (^) is given 
in the Appendix. It shows that with growing N the net 
charge of DNA 77* experiences a first order transition from 
negative to positive values. Here we concentrate only at 
large enough N, where rf is positive. 

Let us make two simplifying approximations. Firstly, 
we assume that the concentration N of PE in the so- 
lution is large enough so that PE translational entropy 
term (the last term in Eq. (0)) can be neglected. In 
other words, we calculate the maximum possible charge 
inversion. This limit is reached when N 3> Nq, where 



Nq — Vq 1 exp(-Zl B \n{C/a)/a) 



(5) 



is an exponentially small characteristic concentration. 
For a long PE Nq is so small that one does not need 
alarge N to get to this limit. 

Secondly, as a good approximation, one can now ne- 
glect second term of the right side of Eq. (||), which is 
responsible for the entropy of defects on DNA. This easily 
leads to a solution for the net charge density 



V 



ln(£/a) 



2a ln[(r s + R)/R] 
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Now one can check that this solution is consistent with 
the assumption that the entropic term can be neglected 
by substituting it back into Eq. (|J). 

Equation (0) shows that 77* is positive indicating that 
the bare DNA charge is inverted. Knowing 77* and using 
1 770 1 = 0.6e/A~ 3.9e/o the charge inversion ratio can be 
calculated 



0.13 
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For DNA R = IOA and a = 6.7A, so that at r s > 10 A 
the ratio of logarithms can be only slightly larger than 
unity. Thus, the charge inversion ratio created by frac- 
tionalization is limited by 20%. Up to such point we 
indeed can neglect Coulomb interactions between defects 
in the chemical potential of the gas of defects (the first 
term in the right hand side of Eq. (|l|)). 



Remarkably, the extremely crude bead-and-stick 
model of PE discussed above can give reliable and uni- 
versal predictions. The calculation described above is 
not sensitive to many microscopic details and chemically- 
specific effects on atomic scale. One could worry about 
behavior of dielectric constant of water at small dis- 
tances, destruction of water solvation shells, other in- 
teractions (van der Waals, hydrogen bonds, etc.) All of 
them are not important because they all modify energy 
of interaction of PE with DNA which does not enter in 
the above calculation. This energy is identical for con- 
figurations on Fig. [l]a and Fig. [ljb. The only difference 
between these configurations is the self-energy of a free 
PE molecule, which does not depend on any details of the 
PE-DNA interaction. Only this self-energy drives charge 
inversion. 

One could also ask about the role of the finite flexibility 
of PE for the tails. As we all know, freely jointed chain 
model of polycation is useful on length scales of several 
nanometers, but is not literally valid even on length scales 
of 6-7 A. We want to emphasize that we do not need ideal 
flexibility of tails, which lets them to be perpendicular to 
DNA cylinder surface. The only requirements for flexi- 
bility of tails assumed in our calculation is that the tail 
can be raised in such a way that its end monomer avoids 
the end monomer the neighboring PE molecule. This 
requirementSpis fulfilled in many cases, for example, for 
the sperminetLS (Z = 4). (One should take into account 
that the neighboring PE charged monomers are usually 
connected by a chain of several neutral monomers). 

Small arches shown on Fig. |l|b, however, are more 
sensitive to flexibility than tails. If the persistent length 
of PE, I is larger than the distance between charges, a, 
loops (arches) have a typical length I. In a long PE where 
arches dominate this leads to replacement of ln(£/a) by 
hx(C/l) in Eq. (Q) and therefore to a somewhat weaker 
charge inversion. 

In the mostly theoretical case of a short and extremely 
rigid PE when even tails can not bend at all, so that a 
PE charge of a neutralized DNA is totally incompress- 
ible, both fractionalization and charge inversion disap- 
pear. This is similar to what happens when Z-ions are 
hard spheres and one layer of thera-exactly compensates 
the uniformly charged backgroundE3. Charge inversion 
disappeares in these cases, because there are no inter- 
nal degrees of freedom of molecules to make the system 
compressible. 

Until now we talked about one-dimensional periodic 
chain of negative charges. If we recall that in DNA this 
chain actually is a spiral we face another requirement 
for the flexibility of a long PE. A PE molecule should 
be flexible enough to follow DNA spiral. Most of PE 
can do that, for example spermine doesE3. On the other 
hand, extremely rigid long PE can not follow a spiral 
of charge and, therefore, screens DNA as an uniformly 
charged cylinder, namely PE rods in this case arrange 
themselves at its surface collinearly with the cylinder axis 
and each other. 
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FIG. 2: The origin of charge fractionalization for a PE with 
linear charge density twice larger than for a DNA spiral, a) 
A Wigner-crystal-like ground state of a periodic chain of neg- 
ative charges neutralized by PE molecules with Z = 3 and 
b — a/2 (their charges are shown by black spheres). A new 
PE molecule is approaching DNA. b) The new PE molecule 
is "digested" by DNA. Its charge is split in +e charges of Z 
grain boundaries. 
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FIG. 3: The origin of charge fractionalization for multivalent 
counterions. a) A Wigner-crystal-like ground state of a peri- 
odic chain of negative charges (white spheres) neutralized by 
multivalent counterions with Z — 3 (larger black spheres) . b) 
The new Z-ion is "digested" by DNA. Its charge is split in +e 
charges of Z grain boundaries. Broken circles are explained 
in the text 



Concluding this section, we would like to say that the 
discreteness of charges does not prevent charge inversion 
even in the worst case of perfect matching. 



III. POLYELECTROLYTE WITH 
NON-MATCHING DENSITY OF CHARGE 

How does fractionalization work when distance be- 
tween charges of PE, b, is not equal to the distance be- 
tween charges of an unfolded DNA spiral, a? Consider, 
for example, commensurate PE with b = a/2, which has 
linear density of charge twice larger than a DNA spiral. 
In this case, PE molecules due to Coulomb repulsion form 
an analog of Wigner crystal where PE molecules alternate 
with vacant places (see Fig. ||a). Even if the PE is ab- 
solutely rigid a new PE molecule creates Z distant grain 
boundaries (domain walls), where one vacancy is missing 
(see Fig. ||b). The charge of each grain boundary is +e, 
so that charge of the new PE molecule is fractionalized. 
and a part of the self-energy of PE is eliminated in the 
way similar to what happens in the case of matching PE. 

Fractionalization continues to work when the linear 
charge density of multivalent counterion (Z-ion) is even 
larger. We can imagine such limit, when replacing PE 
with a metallic multivalent ion (for example, La +3 ), 
which touches only one negative charge of DNA. Then 
we arrive at a ground state of neutralized DNA which 
resembles Wigner crystal even closer (see Fig. ||a). Frac- 
tionalization of a new charge into Z monovalent charges 
of grain boundaries (see Fig. ^b) decreases self-energy 
and drives charge inversion. In contrary to obvious elim- 
ination of interaction between monomers in the case of 
Fig. |l} it is more difficult to see how self-energy is elim- 
inated in Figs. |^ and |^, For example, to get an idea how 
this happens in the case of Fig. || it helps to draw a 
sphere with radius a bit larger than a both around the 



new free Z-ion on Fig. |3ja and around the center of each 
of Z domain walls on Fig. |^b. They are shown by broken 
circles. Let us consider now what happened to the energy 
of the electric field of the new Z-ion concentrated in the 
external space of these spheres. Due to fractionalization 
of Z-ion the energy of the charge Z is clearly replaced by 
the smaller sum of Z energies of monovalent ions. This 
illustrates what we mean talking about elimination of the 
self-energy in this case. 

In Fig. H we already arrived at a model of charge 
density wave and fractionalization in polyacctilcnc and a 
very crude picture for the fractional Hall effect at filling 
factor 1/3. In the latter case, empty circles mean discrete 
Landau states and an electron charge — e is split in 3 
charges — e/3. 

Fig. U also resembles what happens in the case of ad- 
sorption of Z-ions on the line or.-surface with uniform 
distribution of background chargeoli3. In that case, the 
charge of a new Z-ion is smeared along the background 
due to small elastic deformations of Wigner-crystal-like 
strongly correlated liquid. In other words, Z-ion is frac- 
tionalized into infinitesimally small portions. One can 
visualize the transition to the case of uniform surface 
charge imagining that both elementary charge of our lat- 
tice and lattice constant a vanish, while charge density 
of DNA and charge of Z-ion are kept constant. 

Let us return to adsorption of PE with a finite linear 
charge density on DNA and discuss more complicated 
situations, when b < a, but b and a are incommensurable. 
Even in this case ground state of a neutralized DNA is 
a crystal. If an additional PE molecule is adsorbed it 
is still fractionalized to Z grain boundaries with charge 
+e. The only difference from commensurable case shown 
on Fig. U is that grain boundary can include several PE 
molecules. 




FIG. 4: The origin of a shear deformation string between two 
+2e defects, which appeared as a results of adsorption of an 
additional Z-ion with Z — 4 in a two-dimensional generaliza- 
tion of the model of Fig. ^[ 



IV. FRACTION ALIZATION OF 
POLYELECTROLYTE CHARGE IN 
TWO-DIMENSIONS 
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FIG. 5: Fractionalization of an additional PE molecule with 
Z — 3 into three positive defects at a two-dimensional lattice 
of negative charges (white spheres) neutralized by already ad- 
sorbed PE. Positive charges of PE molecules are shown by 
black spheres. 



It is well known that fractionalization of charge of Z- 
ion into free (to move to infinity) grain boundaries shown 
for example in Fig. [}] can not be generalized to a two- 
dimensional case. Let us imagine a two-dimensional ana- 
log of the problem of Fig. a using a square lattice of 
monovalent negative charges which is neutralized by Z- 
ions with charge Z = 4 forming a square lattice with 
period 2a. If we bring another Z-ion and try to split 
it into two point like grain boundaries with charge +2e 
along the main axes of the square lattice, we realize that 
the lattice of Z-ions looses energy everywhere between 
them because of the shear deformation created (See Fig. 
[1J). This is equivalent to a string between charges +2e 
with energy proportional to length. Thus, charges +2e 
of the defects are not permitted to move very far away 
from each other. In other words, they are confined in a 
finite domain. 

Nevertheless, in the first approximation, charge inver- 
sion can be still calculated as if products of the Z-ion frac- 
tionalization were free to move to infinity. Indeed, one 
can estimate the defect confinement size and find that it 
is much larger than a\/~Z (the average distance between 
Z-ions on the surface), because the energy of abovemen- 
tioned strings is proportional to the shear modulus of the 
Coulomb lattice of Z-ions on the negative background 
lattice which is known to be numerically small. There- 
fore, most of the self-energy of Z-ion concentrated in the 
electric field at radius larger than a\[Z is eliminated in 
spite of defects confinement. 

Remarkably, for a reasonably flexible matching PE 
fractionalization into free tails and arches is not a strictly 
one-dimensional phenomenon. It is easy to see that 
the same mechanism applies equally well to a two- 
dimensional square lattice of discrete negative charges 
with the lattice constant, a, equal to the PE bond length 
b. Indeed, one can see in Fig. | that all previous argu- 



ments about the role of tails and arches can be carried 
over to this case. There are no strings between tails and 
arches in this case. This is a remarkable consequence 
of the involvement of additional degrees of freedom re- 
lated to the third dimension. We do not know any other 
classical example of charge fractionalization in a really 
two-dimensional system. 

There are, however, small modifications of the ana- 
lytic formulae for charge inversion. Defects with +e 
charges form now a two-dimensional gas with concen- 
tration cr*/e, where a* is the net positive surface charge 
density playing the roles of rf. The chemical potential 
of this gas is fcsTln(a 2 (T*/e). The surface potential is 
i/)(0) = 2ira*r s / D. The balance of the chemical poten- 
tial of PE molecules adsorbed at the surface with that of 
a free PE in the solution reads 

^4^ + ^ln4, + ^l«. (8) 
D Da a e a a Ze 

Again, assuming that the PE concentration N is large 
(or calculating the maximum possible charge inversion) 
the solution to Eq. (||), for a ~ Ib, within a numerical 
factor, is 

a* ~ (e/or.)/ln(r,/o) . (9) 

One can see that, for a ~ Is, in the free energy gained by 
fractionalization of the PE molecule charge, the entropy 
contribution is comparable to the self-energy, in contrary 
with the one-dimensional case, where the entropic term 
can be neglected. This is due to a higher number of 
degrees of freedom which a two-dimensional surface pro- 
vides to the gas of defects. If r s ^> a, the charge inver- 
sion ratio for the two-dimensional case is smaller than for 
DNA: 

= — In — . (10) 

e/a z r s a 
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An important role of elimination of the self-energy for 
adsorption of a flexible PE on an oppositely uniformly 
charged surface can be traced in Refs. 0, 0, [ll| 



V. CHARGE INVERSION OF DNA BY 
DENDRIMERS AND FRACTIONALIZATION 

Until now we considered adsorption of linear charged 
molecules (PE) both on one- and two-dimensional lat- 
tices of the background charge. It is interesting to note 
that the fractionalization mechanism works for molecules 
of other shapes, too. Let us, for example, consider den- 
drimers (star-like branching molecules with a large num- 
ber of monovalent positive charges on their periphery)., 
which were also shown to invert the charge of DNAtS. 
Dendrimers with charges ^=4, 8 can easily compensate a 
compact group of nearest Z charges of both DNA helices. 
If a DNA double helix is totally covered and neutralized 
by such dendrimers (see the schematic Fig. ||a for Z=A) 
an additional dendrimer can still be adsorbed on DNA. 
This happens because two charges +e of two distant al- 
ready adsorbed dendrimers can be raised above the DNA 
surface when a new dendrimer molecule is adsorbed on 
DNA (see Fig. ^o). As in the case of linear molecules, 
this fractionalization of the dendrimer with charge +4e 
into two charges +2e leads to the gain of its self-energy 
and to charge inversion. 

Again we see that all these phenomena became possible 
only due to the additional of freedom of PE molecules, 
which in this case is rotational. (Fractionalization into 
charges +e in this case can leads to a larger energy be- 
cause all adsorbed dendrimers between two dendrimers 
raising one tail should be deformed leading to a string 
with energy proportional to the length between them.) 

If we deal with higher generations of dendrimers which 
have very large charges such as 32e or 64e, we arrive 
at a different Wigner-crystal-like picture (see Fig. [?]). 
Because of the three-dimensional structure of their chem- 
ical bonds these molecules can not expand enough so that 
each charge of them reaches an opposite charge of DNA 
and compensates it. In other words, when projected to 
a DNA double helix, these high generation dendrimers 
have much larger linear density of charge than the double 
helix itself. Thus, large segments of the helix between ad- 
sorbed dendrimers remain negatively charged, and form 
a Wigner-Seitz cells around each dendrimer. This is 
how with growing charge of dendrimers the fractionaliza- 
tion mechanism is replaced by the mechanism of Wigner- 
crystal-like correlations. Qualitative difference between 
DNA complexes with dendrimers of low and high gener- 
ations has been clearly demonstrated experimentallyEj. 
Because large fraction of DNA charges is not neutralized 
by dendrimers the high generation complexes are more 
sensitive to the salt concentration. 




FIG. 6: The origin of charge fractionalization in dendrimer 
adsorption, a) Two linear chains of negative charges of DNA 
(white spheres) which are obtained unfolding DNA spirals of 
charges. They are completely neutralized by positive den- 
drimer molecules with Z = 4. Dendrimer are schematically 
shown by plane crosses with +e charges (black spheres) at the 
ends. A new dendrimer molecule is approaching DNA. b) The 
new dendrimer molecule is "digested" by DNA. Its charge is 
split in +2e charges of the tail doublets. 



VI. SHARING OF DNA CHARGES AS A 
MECHANISM OF CHARGE INVERSION 

Let us return to complexation of a DNA double helix 
with PE molecules with the matching bond length, b = a, 
and discuss another possible mechanism of charge inver- 
sion, which is also related to the discreteness of DNA 
charge and further increases the positive charge of DNA- 
PE complex. Let us consider a monomer tail of PE on 
Fig. [l]b and explore whether some energy can be gained 
if the positive charge of this monomer moves down to 
the plane of DNA charges, approaches already neutral- 
ized negative charge of the DNA and shares it with the 
end monomer of the neighboring PE molecule in a way 
shown in Fig. |^. If these two end monomers may sit 
on exactly opposite sides of the negative charge of DNA, 
the additional energy e 2 /2c? can be gained, where d is the 
distance of the closest approach of a PE monomer and 
a DNA charge. At a sufficiently small d this energy can 
be even larger than the gain per tail from elimination 
of the self-energy. In a DNA double helix, all the neg- 
ative charges indeed are on the ridge above neighboring 
neutral atoms. Two sufficiently small monomers may fit 
into the large and small groves on both sides of the ridge. 
On the other hand, if because of sterical limitations they 
can not be in the perfect opposition the energy gain is 
smaller. If both end monomers PE have the same size 
as the negative charge of DNA the additional energy of 
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cussed in a model of uniformly charged macroion surface, 
but still is accompanied by fractionalization. 

In conclusion, we emphasize that in any case charge in- 
version happens due to the fact that a new PE molecule 
rearranges already adsorbed PE in such away that its 
image or correlation hole strongly attracts this new PE 
molecule. This physics can not be described by the 
Poisson-Boltzmann theory because this theory uses the 
mean-field potential which does not depend on the posi- 
tion of a new PE molecule. 
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FIG. 7: The origin of charge inversion in adsorption of high 
generation dendrimers. a) DNA double helix (gray) neutral- 
ized by a Wigner-crystal-like liquid of a high generation den- 
drimers (dark spheres). A new dendrimer molecule is ap- 
proaching DNA. b) The new dendrimer molecule is integrated 
into Wigner-crystal-like liquid while neighboring already ad- 
sorbed dendrimers slide away from it and smear its charge 
over the helix. 




FIG. 8: A view from the top on an unfolded spiral of nega- 
tive charges of DNA (white spheres) and two PE molecules 
(black). Two positive end monomers share a negative charge 
of DNA in the perfect opposition. 



sharing vanishes when all three spheres touch each other 
forming equilateral triangle. This still leaves room for 
sharing effect, while say one monomer perfectly fits in 
large grove but the second one only partially fits in the 
small grove. 



VII. CONCLUSION 

The main result of this paper is that discreteness of 
surface charges of a macroion, for example, double he- 
lix DNA does not prevent its overcharging by oppositely 
charged PE and other Z-ions. For a flexible PE even 
in the worst scenario of matching PE and DNA geome- 
tries, when in the neutral state all charges of DNA are 
perfectly neutralized by PE, charge inversion happens 
due to fractionalization of PE charge into +e charges 
of defects. This is extremely transparent mechanism as 
illustrated on Fig. [l| It is clearly related to internal de- 
grees of freedom of a flexible PE. In the non-matching 
cases, the mechanism of charge inversion for discrete sur- 
face charges looks more similar to the one previously dis- 
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APPENDIX A: CHARGE OF DNA AS A 
FUNCTION OF POLYELECTROLYTE 
CONCENTRATION 

In Sec. 0, we assumed the bulk PE concentration, 
N, is very large so that the translation entropy cost of 
condensing them on DNA can be neglected. Equation 
@ , thus, gives the upper limit for the DNA net inverted 
charge. On the other hand, at very small N, the entropy 
cost cannot be neglected and leads to the undercharging 
of DNA. In this appendix, we would like to calculate 
rf (N) explicitly and show that in the matching case DNA 
molecules change their sign with increasing N by a first 
order phase transition. 

When DNA is undercharged, rf < 0, instead of a gas 
of raised monomers (tails and arches) on the DNA sur- 
face, one has a gas of vacancies. These are the DNA 
charges which are not covered by any PE monomers. At 
low concentration (small undercharging) , these vacancies 
practically do not interact and their chemical potential 
can be approximated by that of an ideal gas at the same 
concentration fcsTln(?7o/|?7*|). Thus, Eq. (jij) needs only 
a small modification to properly describe both the over- 
and under-charged DNA: 



2,* ^ 



e/l 



+ R . m 

m -, r 

R \v*\ 



Z N 



(Al) 



where the second term is the chemical potential of raised 
monomers in the overcharging case and it is the chemical 
pote ntial of vacancies in the undercharging case. In Eq. 
( Al), we have also combined two terms of Eq. (Q) using 
the characteristic concentration No given by Eq. (||). 

It should be noted that, the apparent divergence of the 
left side of Eq. ( |Al| ) at small rj* is related to the fact that 
we neglected a small concentration of intrinsic defects 
(raised monomers and vacancies). This concentration is 
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a) 



b) 




N C N, 



ln(N/N )/ Z 



\n(N/N Q )/ Z 



FIG. 9: Behavior of rj*(N) near the charge inversion transi- 
tion. The limit discussed in Sec. [n| starts only in the u pper - 
right corner of the figures, a) Graphical solution to E q . (Af ). 
The solid line is the right hand side of Eq. ( |Al| ). When 
N > Ni, there is only one solution 77*. When N < Ni, there 
are three solutions, 771,2,3- b) The DNA net charge density if 
as a function of N (the solid line). A first order phase tran- 
sition from the undercharged to overcharged state happens 
at N = N c . The dashed lines correspond to the metastable 
values of rf . 



of the order a -1 exp(— e 2 /2dksT). It exists even at 77* = 
and truncates the divergence of ln(ryo/|??*|) at small 77*. 
However, when \rf\ 3> ea' 1 exp(— e 2 /2c?fcsT) one type of 
defects dominates over the other and one can neglect the 



contrib utio n from the minority ones. This is what we did 
in Eq. (|Al]). 

To understand how 77* varies with N, it is very in- 



structive to solve Eq. (Al) graphically. One can see the 
following behavior: 

When N is large such that N > Ni, where 
i/N )/Z = -1 - ln{2ijoln[(r. + R)/R]/(e/l B )}, Eq. 
has only one solution for 77*. This solution is posi- 
tive, indicating that the DNA helix is overcharged. 

When N decreases a little bit belowJVi , there are three 
solutions 77* > > 77* > 77* of Eq. Q (Fig. |a). The 
solution r/l corresponds to the stable overcharged state. 
The solution 773 corresponds to the metastable under- 
charged state. The solution 772 is unstable (it corresponds 
to a local maximum in the grand potential of the system 
located between two local minimums at r/l and 773.) 

When N decreases below N c where N c is defined as the 
PE concentration at which the two shaded areas in Fig. 
|9|a equal each other (Maxwell rule), a first order phase 
transition happens. (Note that at exp(— e 2 /2dfcsT) <C 1 



a calculation of N c can be done with help of Eq. (Al) 



because the truncation due to intrinsic defects produces 
only a small correction to one of areas.) The overcharged 
solution r/l becomes metastable while the undercharged 
solution 773 becomes stable. Thus, the function 77* [N) 
has a finite jump at N = N c . This function is plotted 
by the solid line in Fig. [|b. The metastable branches of 
77* (iV) are plotted as the dashed line. 
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